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We study the stress energy two-point function to show how short distance correlations across the 
horizon transform into correlations among asymptotic states, for the Unruh effect, and for black 
hole radiation. In the first case the transition is caused by the coupling to accelerated systems. 
In the second, the transition is more elusive and due to the change of the geometry from the near 
horizon region to the asymptotic one. The gradual transition is appropriately described by using 
affine coordinates. We relate this to the covariant regularization used to evaluate the mean value 
of the stress energy. We apply these considerations to analogue black holes, i.e. dispersive theories. 
On one hand, the preferred rest frame gives further insight about the transition, and on the other 
' hand, the dispersion tames the singular behavior found on the horizon in relativistic theories. 
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I. INTRODUCTION 

O ' In Minkowski space the vacuum of quantum fields, the ground state of the Hamiltonian, is globally defined. For 
instance, for a massless free field (j) in two dimensions, it is the state annihilated by the destruction operators 

^ \ associated with the positive norm modes 0fe = e~*^*+'^^/(47ri7)^/^ where = \k\ > 0. When probed locally, it 
determines correlation functions, such as that of the stress-energy tensor. Considering Tuu = {du(t>Y, where U {V) 
is the retarded (advanced) null coordinate U ^ t — z, {V = t + z), the connected function is 

^ . {Tuu{U)Tuu{Uo)) - (duduo{4'{U)^{Uo))] 
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[47^{U-Uo-^e)^) ' 

, The regulator e — )■ 0+ specifies the nature of the singularity in the limit U ^ Uq. It arises from the fact that only 
■ positive frequency ft = idjj contribute to Eq. ([ij. 

Even though this function monotonically behaves in AU^^, it determines subtle effects. This can be seen by 
considering accelerated systems where both local and non-local phenomena related to the Unruh effect [11-01 can be 
. ^ obtained from Eq. ([iJ . Moreover the correlations of Tjji; in the (Unruh) vacuum across the future horizon of a black 
, hole behave as in Eq. ([T]). However since the black hole geometry is non uniform, as one probes the vacuum further 
^ ' away from the horizon, there is agradual change which results, on one hand, in correlations typical of a thermal fiux 
5^ [ at the Hawking temperature 0, U |^ and, on the other, in correlations across the horizon between Hawking quanta 
and their partners [Ifll ■ 

The aim of the present paper is to analyze the black hole case in the light of the Minkowski case. We show that 
the gradual change in the correlations comes through the observables used to probe the state, and not from 4> which 
is conformally invariant and thus insensitive to the geometry. We relate this gradual change to that obtained when 
considering renormalizedj;xpectation values which also depend on the conformal factor. In the last part, we consider 
dispersive field theories [ll!], or analogue black holes in condensed matter [T^. [isj. First, we show that the preferred 
rest frame [l^ associated with dispersion defines new scalars (observables) which give an unambiguous meaning to 
the gradual change of the stress energy correlations. Second we study how dispersion affects the correlation pattern.^ 
In addition to the robustness of the spectral properties of Hawking radiation [ll], HI, fl^- that of the long distance 
correlations across the horizon was established by constructing wave packets of in modes [Ig] . We shall return to this 
observation and make contact with the analysis of density-density correlations in Bose- Einstein condensates [18, - .20j . 



* Electronic address: renaud. parentani@th.u-psud.fr| 

^ While completing this work, we became aware of [15| where similar conclusions are reached. 
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We have organized the paper as fohows. In Section II, we describe the behavior of the stress energy 2pt function 
in flat space, and interpret its properties by introducing external (accelerated) systems which probe the field. In 
Section III, we perform the same analysis in a black hole geometry. Using affine coordinates, we present an invariant 
description of the transition from its near horizon behavior, which is similar to Eq. (jlj, to its long distance behavior. 
Finally, in Section IV, we consider the stress energy correlations in dispersive theories, and investigate what are the 
modifications induced by dispersion. We conclude with remarks about the similarities between these dispersive effects 
and those obtained with relativitic fields propagating in fluctuating metrics [2l|, . 



II. UNRUH EFFECT AND CORRELATIONS ACROSS A RINDLER HORIZON 



We introduce the null coordinates au ^ — ln(— aC/), av = \n(aV), defined respectively for t7 < and V > 0. These 
are related by r = {u + v)/2 to r, the proper time of an accelerated system which follows — ~ I/a^ in R, 
the right quadrant z > \t\. (a is the constant acceleration.) When U and Uq are negative, i.e. with both points 
on the right oi U — 0, the future horizon of the accelerated system, Eq. ([T]) can be re-written using u. Introducing 
= {du4>? = Tuu Eq. © becomes 

(T-(.)r-(uo)) = (2) 

Since l/(sinhy)^ = S„ (y + imT)~'^, where the sum over n goes from —oo to oo, (T'^T^J is the periodic function in 
Imu of period 27r/a whose zeroth term, n = 0, behaves as in the vacuum, see Eq. ([T]). Hence Eq. ^ is the 2-point 
function in a thermal bath at a temperature equal to a/27r. Noticing that along the accelerated trajectory, one has 
T = u (because v = u) Eq. ^ determines the correlation function of the energy density probed by the accelerated 
system. In this way we recover the Unruh effect: accelerated systems perceive the vacuum as a thermal ensemble at 
that temperature 

It is less usual but equally interesting to consider Eq. ([T]) with one point on either side of the horizon. Using u and 
Tuu on the right, and their symmetrical counterpart on the left, au = \n{aU) and Tuu = {du4')^ = Tuu for t/ > 0, 
Eq. ([1]) becomes 

{TL{u)Tt{u,)) = + • (3) 

Quite surprizingly (r^fi^T,^) possesses a 'bump' which is centered at u + mq = 0. In terms of the affine coordinate 
[/, it is located at IJ = —Uq, on the opposite value from the horizon. This is unnexpected because, at fixed AV, 
{(j){U, V)(j){Uo, Vo)) — ^ In AJ7 + Cst. has no maximum when expressed in terms of u and uq. In fact, the maximum 
in Eq. ([3]) results from the combined effect of the Jacobians dU/du associated with the tensorial character of Tuu, and 
the monotonic power law decay of Eq. ([I}.^ These properties apply to massless fields in any space-time dimension 
d. Because of the dimensionality of the cj) field (= {d — 2)/2) the location of the maximum is shifted a.t u + uq = 
— -arctanh^^. 

So far however Eq. ([3]) is simply a re-expression of Eq. ([T}: no physics is gained when writing a tensor in a new 
coordinate system. So, isn't the maximum of Eq. ([3]) just a coordinate artefact ? No, Eq. ([3]), as Eq. ([2]), governs 
physical effects. For instance, when looking at the state of two opposite accelerated systems [7|, the correlation matrix 
of this bi-partite system exhibits a non-trivial entanglement near the maximum of Eq. ([3]) . Similarly the scattering by 
two opposite accelerated mirrors destructively interferes near that maximum [6)] . As another example , expressions 
closely related to Eq. (jSj obtain when looking at the value of which is correlated to the fact that an accelerated 
detector in R has made a transition. In all cases the physics is the same: measurements performed in R are statistically 
correlated to those performed on the other side of the horizon. 

To convey the idea that the re-writing of Eq. ((IJ as Eqs. (HJ [3]) is not merely a coordinate change but has to do 
with the physics of accelerated systems, we make two observations. First it is conceptually useful to introduce the 
{U part of the) scalar energy density that such system in R will measure: = Tf^^u'^u'^, where = dx^/dr is the 
unit vector field tangeant to its trajectory. Then the l.h.s of Eq. ^ is {p^{t) p^(ro)). Similarly the l.h.s of Eq. is 
{p^{f ) p^-{to)). Hence both Eq. ([2]) and Eq. ^ have a clear interpretation as energy-energy correlation functions of 



My attention was drawn to this point in a discussion with W. Unruh and T. Jacobson that took place during the workshop "Towards 
the observation of Hawking radiation in condensed matter systems" held at IFIC in Valencia in Feb. 2009. 
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accelerated systems. Secondly, we notice that Eq. ([2]) only depends on u — uq. This results from the invariance of the 
vacuum under Lorentz transformations u ^ u + b,U — >■ Ue~°'^. In fact the Minkowski distance between two points 
situated along the accelerated trajectory v = u is 

9 A 9 A 9 4 sinh^ (a(T - To)/2) 

a"' 

We therefore see that the stationarity of Eq. in u for arbitrary values of v are related to the properties of s^(At) 
evaluated along v = u. In a similar fashion, Eq. ([3]) is related to the distance s^p bewteen two points situated along 
the opposite accelerated trajectories defined by z"^ — P = Using the proper time tq (f) to localize the point on 

the right (left) trajectory, one gets 

2 _ 4cosh'(a(f + ro)/2) 

The maximum of correlation in Eq. ^ thus coincides with the minimum (space-like) distance bewteen the two 
trajectories. Because of Lorentz invariance no preferred spatial distance Az^ is selected though. 

Yet, in preparation for the black hole problem, one would like to know what is the origin of the mathematical 
properties of Eq. © and Eq. ([3]). In fact, these must be deeply rooted since Eq. ([3]) is the analytical continuation 
of Eq. (21) obtained by subtracting half of the imaginary period that fixes the Unruh temperature. In addition, we 
wish to identify what could explain the maximum of correlation in Eq. ([3]), which is independent of A.V, and thus 
not necessarily refers to Eq. ([5]) and accelerated systems in L. As shown in the Appendices, these properties can be 
understood when using the Unruh modes [l[ . In brief, the maximum of Eq. ([3]) finds its origin in the entanglement 
in Fock space between states with opposite Killing frequency w = idu- This entanglement translates into space-time 
correlations because, in each pair, the partner mode lives on the other side of the Rindler horizon, and is weighted in 
Eq. (jA2p by a factor z^^ which is real for all w. In addition the maximum value of correlations is directly related to 
the mean occupation number of Rindler quanta. 



III. HAWKING EFFECT AND CORRELATIONS ACROSS THE HORIZON 

We have a double aim. First, we wish to interpret the equivalent of Eq. ([T]) and Eqs. (|2l3p in a black hole geometry. 
Second, we aim to analyze the differences between the black hole case and Minkowski. In particular we are looking 
for a smooth and coordinate invariant interpolation from the near horizon region, where Eq. ([l} should make sense, 
to the asymptotic regions where Eqs. (|2I3I) should do. 

To reach these ends, we need to identify what corresponds to the Minkowski vacuum, the affine coordinate U, and 
the u coordinate associated with boosts u ^ u + b,U U e""^. For simplicity we shall work with a two dimensional 
conformally invariant field (f>. Then the identification is rather easy and known [l|, Q . The novelty mainly consists in 
the attention paid to the stress-energy correlation function. 



A. Prelude: conformal invariance 



We work with a conformally invariant field because the simplicity of the expressions will oblige us to identify under 
which conditions identities valid in all space-times acquire a physical meaning related to the Unruh and/or to the 
Hawking effect. 

In two dimensions the line element can always be written in double null coordinates as 

ds^ = ~C{U, V) dUdV . (6) 

In these coordinates, </> obeys djjdvff) — 0- Therefore the conformal factor C drops out and the right and left moving 
sectors remain decoupled = 4>{U) + <P{V). As a result, for any U coordinate, one can consider the {/-vacuum state 
defined by the positive frequency modes (t>n — e~*^'^/(47r51)^/^. In that state, the connected correlation function of 
Tuu = (du^)^ 



{Tuu{U)Tuu{Uo))c = {Tuu{U)Tuu{Uo)) - {Tuu{U)) {Tuu{Uo)), 



(7) 
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obeys Eq. independently of the regularization scheme used to compute {TuuiU)) Introducing the coordinates 
au — \n{—aU) and au = \n{aU) for respectively negative and positive values of U ^ Eq. © and Eq. ^ also follow, as 
mathematical identities valid in all space-times and for any coordinate U . Some physical input is therefore needed to 
transform these identities into meaningful relations among observables. 

To relate Eq. ([2]) to the Unruh effect is straightforward. It suffices to consider a particle detector following the orbit 
u = V. It will perceive the {/-vacuum as a thermal state at a temperature ajlix. (In general a is no longer the proper 
acceleration, consider e.g. de Sitter space where the thermal bath is also perceived by inertial detectors.) Similarly, 
for Eq. when introducing another detector which follows the mirror trajectory u = the combined state of these 
two detectors will be entangled as that of accelerated detectors in Minkowski Q , again as a direct consequence of the 
conformal invariance of </>. This generalization of the Unruh effect applies to any field and in any dimension, but only 
approximatively provided the acceleration is much higher that the space-time curvature. (This limit is also used in 
the thermodynamic analysis of space-time [23j . However the fact that Eq. ([3]) always applies indicates that the purity 
of the quantum state, and not only the mean fluxes, is a key ingredient.) 

To relate Eq. ([2]) and Eq. ([3]) to the Hawking effect is more elaborate. As we shall see, the metric plays a crucial 
role in "selecting" the [/-vacuum and the coordinate u. This follows from the fact that, unlike the connected 2pt 
function of Eq. ([7]), other observables are not conformally invariant, as for instance the renormalized expectation value 
of Tuu Q] 

(W°"--y|^CV2 5^^. (8) 

In Sec. III.D, we shall see that the breakdown of the conformal invariance comes from the renormalization scheme 
and not from the "bare" operator. 



B. The correspondence 

We work with one dimensional stationary black hole metrics that we describe using Eddington-Finkelstein (EE) 
coordinates v, r in which the line element reads 

ds^ = ~{l~w'^{r))dv^ + 2dvdr. (9) 

We shall not consider specific functions for w^, as e.g. 2M/r which describes a Schwarzschild black hole. Hence the 
coordinate r should not be conceived as associated with spherical symmetry. In fact, from a 2D point of view, r in 
Eq. ^ is an affine coordinate along v — est. We also notice that r is affine along r = est., where t is the proper time 
appearing in Painleve-GuUstrand (PG) coordinates [l^ where ds"^ = —dr^ + {dr — wdr)^ . This guarantees that the 
forthcoming analysis can be made in PG coordinates. 

When w is a constant, the geometry is Minkowski. Therefore the differences with Sec. II will all stem from the 
gradient of w{r). We have adopted these coordinates because they are well behaved in the late time portion of the 
space- time which is relevant for the Hawking effect and the correlations across the horizon, see e.g. Q. In particular 
they stay regular when crosses 1, say at r = r^- This locus is the future horizon (for the static observers at fixed 
r > Th) and it is a null line U = est. The absence of the corresponding V (past) horizon induces a disymmetry: left 
moving modes are regular across the horizon U = est. and play no role in the following. 

Assuming that w{r) — >■ Woo = est. for r — >■ oo, the asymptotic observers at fixed r are inertial, at rest with the black 
hole, and their proper time is Tas — v{l — w^)^/^. By rescaling v and r, one can work in a gauge in which Woo = 0. 
From now on we adopt it. The positive norm modes that these observers will use are (jy^ = e' -^""/(47ra;)i/2, with 
w > 0. When r oo, u,v are related to Tas by Tas = {u + v)/2. At finite r, one has u = v — 2 J dr/{l — ■up'). One 
sees that u might diverge when vo^ crosses 1. 

To settle this, one needs to know how w behaves across r^. We assume a regular behavior vP{r) ~ 1 — 2K{r — Th) 
as this is the case for regular collapses, (k > coincides with the surface gravity since Wao = [3(|. The extreme case 
K — requires a special treatment [36j.) At fixed v, one finds u ^ — Kln(r — r^). Therefore, when r — > r'l^^ at fixed v, 
(j)'^ behaves as 

•p-w (4^^)1/2 ' y'^) 



^ The connected 2pt function of T^i/ governs gravitational back-reaction effects beyond those included in the semi-classical Einstein 
equations. [23| This provides its dynamical relevance. 
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and can be taken to vanish for r < rh- The function ry{u) obeys 

dry — ——C{u,v)du, (11) 

where C{u,v) = (1 — w^) is the conformal factor in the u^v coordinates. It should be emphasized that ry{u) is an affine 
parameter for aU w{r), and this property is governed by C which appears here for the second time. (It was imphcitly 
used before when imposing that the coordinates u, v are related to the asymptotic proper time by 2Tas ^ u + v.) 

When replacing the acceleration a by k in Eq. (jlOp . the correspondance with Eq. (jA2p is manifest, and physically 
meaningful because r„ is affine, as U is in Minkowski space. However, this correspondence is confined near the horizon 
since (f>f;f behaves as e'^"'' for r — > cx) at fixed v. Therefore, the gradual change of the mode (j)^ from ^ (r — r^)"^/** 
to ~ e'^"'' is what distinguishes the back hole case from Minkowski. These facts have been recognized in P, [^J] and 
shall be further exploited below. 

The second part of the correspondence concerns the state of (p. While in Minkowski space the notion of vacuum is 
unambiguous, in a black hole geometry this is lost. Nevertheless the late time behavior [8| is universal and stationary 
when described in terms of u (or Tas). This can be verified by considering different collapses and different initial 
states: in each case there are transients, but after a few e- folding it-times 1/k, these fade out as ~ e"''" and the 
stationary values set in (unless the collapse and/or the state is singular). Therefore, as far as the description in terms 
of u is concerned, a single stationary state is selected. As noticed in [l|, this state is most simply characterized by 
Unruh modes, (j)uj, exactly as the Minkowski vacuum can also be, see the discussion after Eq. (|A2[) . The particular 
combination of modes weighted by = means that only contains positive frequencies Ox = iduK where 

Uk is the regular coordinate across Th which is related to u by dUxIdu — e"''". Indeed, stationarity requires an 
exponential relation bewteen u and Vk^ whereas regularity across fixes the decay rate dlnUK/du to be k. 

The correspondence is now completed: the regularity and the stationarity of the state allow to characterize it by 
the modes of Eq. (jA2p . with U replaced by Uk = —e~'^"/K. This exponential is no longer related to a Lorentz 
transformation, but is still related to an isometry whose Killing field norm oc (1 — w'^Y^^ vanishes on the horizon. 

C. Correlation functions and asymptotic quanta 

In the Unruh vacuum, the connected 2pt function ofTuu = (Sfy^^c^)^, see Eq. ([7]), is exactly given by Eq. ([Ij. Then, 
since nu = — ln(— K[//f), Eq. ^ with a — > k also obtains. However, as discussed in HI. A, there is nothing special 
about these identities. What makes Eq. ([2]) meaningfuU here is that both Uk and u are affine, respectively across the 
horizon and asymptotically. Hence Eq. ([2]) tells us that when evaluated in the regular state, the 2pt function used by 
inertial asymptotic observers is that of a thermal flux at the Hawking temperature k/2tt |3|. 

Let us now consider the equivalent of Eq. ^ . Introducing on the other side of the horizon ku = In kUk , Eq. ^ with 
a — >■ K automatically obtains. But what does it mean ? Unlike for the Unruh effect, we cannot refer to accelerated 
systems in the L quadrant. Nevertheless the procedure of applies to black holes, and one can study the T^i^ 
correlated to the detection of an asymptotic quantum. As in Minkowski, there is a reduction of the state: expectation 
values should be computed with the reduced density matrix, see App. C. Then, because of the entanglement of 
Eq. (|B2[) . there is a correlation between T^j, evaluated inside the horizon and a detection on 1+ [l3|. Equivalently, 
one can directly look for correlations in energy across the horizon and obtain Eq. 

There is however an important difference between Eq. ^ and Eq. It originates from the different status of the 
coordinates u and u. Eq. ([2]) has a clear interpretation because du at fixed v is affine for r — )> oo. If we can consider 
the equivalent of Eq. Q with both points inside, and obtain the same answer when using u, we should ask under 
which conditions would (T^u'^uu) posses an intrinsic meaning (without any reference to the external region). For this, 
it "suffices" that the inside region be also infinite and that w — > const, for r — — oo. Whereas it is unlikely that this 
be relevant for black holes, we assume it is the case and consider the consequences. (In a next Section we discusss 
analogue black holes where this possibility can easily be realized.) When w — >■ const, for r — >■ — oo, the notion of 
asymptotic quanta equally applies to the negative frequency partners [35| . Then, {T^y^T.^y} ^^^^ the same meaning as 
{TuuTuJi and the long distance correlation {T^yT^y} of Eq. (O can be probed. 

D. Locality, covariance and renormalization 

To interprete Eq. ([7]) we have so far used the asymptotic properties of the metric. However we would like an 
intrinsic description of the gradual transition from the horizon to the asymptotic regions. To this end we consider 
the renormalization procedure. As we shall see its covariance supplies the intrinsic description we are looking for. 
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and this by breaking the conformal invariance. (It should be clear that different rules on how to interpolate wih give 
different behaviors since the 2pt function of Tjjjj is not a scalar) . 

We noticed that both u and Uk are afhne respectively for r — >■ cx) and r ^ rh- We also noticed that r„ of Eq. (fTTj) 
is afhnc all the way through. In what follows we exploit this to relate the 2pt function of Tjju to its renormalized 
value of Eq. (|8]). To this end, we study the coincidence point limit of Eq. (O in the Unruh vacuum using 

Trr - {dAf = Tuu, (12) 

which is "coordinate invariant" since is, up to a scale, globally defined. To present the concepts we first work at 
spatial infinity. There, in the limit Uq — >■ u, one has 



/rr f \rr f ^\ f dUKiu) dUK{uo)Y f 1 
{Tuu{u)Tuu[Uo)) K = — 1 1 T 



1 



\ du duQ J \4iT {Uk{u) — Uxiuo) ~ ie) 



1 V , 2(T„„(m)) 



rcn 
K 



. , ■ ; ■ w • ^9 , (13) 

^47r (u — Uq — le) / 47r(u — uq — ley \u — uq J 

where (r„„(M))5^" is the asymptotic expectation value of r„„ in the Unruh vacuum. Indeed, it is defined as 

(r„„(u))^" = lim 9„a„J ((/>(u)0('Uo))k - (0(w)0(wo))m) 

12^ \ du ) du J ASir' ^ ' 

where {ip{u)(p{uo)) m — ^ ln(u — uq — it) is evaluated in the asymptotic (Boulware) Minkowski-like vacuum. The 
procedure of Eq. ([T3|) is clear: by subtracting the asymptotic vacuum divergence, one extracts the excess in the "noise" 
and thus identifies the mean value of Eq. (jl4[) . Moreover, this procedure can be univocally covariantized and applied 
at every space time point. This is achieved by making use of r„(u): 



{Trr{u)Trr{uo))K = 



/I 1 



\47r {rv{u) - rv{uo) + ie) 



where {Tj.r{u, w))^" is the renormalized value of T^r which is defined as 

(Trr(u,w))5?" = lim 9r9ro( (0(w)(?!)(uo))k - (0(M)'/'(wo))local 

Uq^U \ 



_±_ (dUKV'\2 (dUK\"^ 
12^ \ dr ) ''{drj 



(16) 



where 



((?!)(u)0(mo)) local = -y- ln(r„(u) - r^{uo) + ie). (17) 
47r 

In virtue of the covariance, when using the affine coordinate r.u, the subtraction term possesses this universal form. 
This explains why Eq. generalizes Eq. at every point. 

Several remarks should be made. Firstly, using Eq. ([TT|) . one verifies that Eq. p6)) gives back Eq. (|8]) in all space- 
times and all vacua. Hence Eq. (jl6p can be seen as an alternative expression for it. Secondly, since the vacuum 
is defined through the coordinate Uk, and the subtraction only refers to the metric through Eq. (|lip . {T^rY^^ is 
governed by dUx/dr and nothing else. Thirdly, even though {TjjuTuij) k is independent of w, the covariance of the 
divergent terms in Eq. (ITSl) unambiguously defines the i;-dependence of ('r^™)i<- in Eq. ([Tl 
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E. Covariant description of stress tensor correlations 



Following the same logic, we use Trr of Eq. (|12p to characterize the stress-energy correlations in the Unruh vacuum 
at every point. Using Eq. (IT|), Eq. ([7]) gives 



{Trr{v,r)Trr{vo,ro))K = { drdro{<j)iv , r) <l){vo , ro)) K 



/dUi^dUKyfj_ 1 



(18) 



This bi-tensor field depends on AU]^^ since the state is the Unruh vacuum, but unlike Eq. ([T]), it depends on the 
actual location of the two points through the w-dependence of the Jacobians dUx/dr. Using UK{v,r) — e^'^'"UK{r) 
which follows from the stationarity of the metric, we can extract this w-dependence and obtain 

(r..(«, r) T„{v„ ro))K = (^^) ^ [^^{e^^ Uk - - ^e) . (19) 

There is yet another interesting way to write this correlator. Using 



which follows from Eq. (|lip . we get 



{T„{v,r)T„{v,,r,))K ^ (- ,^ ^^^^ ^- 3) ■ (21) 

\^ TT (1 - w^){l - w^) {jJk -U%- ie) J 

From Eq. (jl9p . it is clear that the correlator is a function of only 3 variables, r, tq and v — vq, the state (Unruh 
vacuum) being stationary. Because it still depends on 3 variables, its behavior in different 2 dimensional sections 
illustrates different aspects of the correlations associated with the Hawking effect. When fixing rg, on I+, Eq. (fT8|) . 
function of r, v, describes the correlations associated with a late detection on 1+ [l^. Instead, at equal time vq = v, 
Eq. (flS)) describes the correlations in the r, rg plan that have been accumulated in the past of that time [l^ . 

Before studying these two cases, we need to be more precise about the black hole geometries we shall work with. 
As explained before, we consider profiles that become constant for r — ±00. To have a simple example at hand, we 
choose directly Uk {r) since it is the only relevant function in Eq. ([TQ]) : 



kUk{x) = - (e^'^^ - e-^-^^) , (22) 

where x = r — r^- This is a kind of symmetrized version of the Schwarzschild metric where Uk = — xe^"^. Using 
Eq. ((201), one gets 

1-^2 = 4 o- ■ (23) 

Near the horizon, for <C 1, one has vo^ — 1 — 2kx for all values of k, and the asymptotic values are Woo = for 
X — > 00, and w't^ = 1 + k/k for x — > —00. 

1. Correlations to a late detection 
We fix xo,vo on I^, on the future right infinity, and label it with uq- We consider Eq. ([T8| 

Trr{r,v)\uo.I+ = {Trr{v,r)Trr{uo,Ipj)K, (24) 

as a Ipt function. As such Trr(r,v) is a special case of the conditional value 

frr{r,v)\na = {Trr{v,r)Ilj-+)K, (25) 
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Figure 1: The stress- energy conditional to a late detection. We represent Trr of Eq. (|24p in the x,v plane {v being vertical) for 
a detection at xq = 4, = 0, outside the picture in the top right-hand region where w(r) is constant. The horizon is at a; = 
and K = K = 1 in Eq. (|23|l . When \x\ > 1, on both sides, the pattern is translation invariant along the null direction because 
w is constant. Instead for |x| < 1, one sees the endless focusing of the null lines for u — >■ — oo. The "post-selected" partner 
propagates along v + 2x = —8 = vo — 2xo, i.e. along the opposite trajectory fixed by xo,vo. These features were found in [l^ 
considering Eq. (|25p for a typical Hawking quantum. 

where 112-+ ^ projector that specifies a state on I^, see App. C for a brief account, and App. C of "s*] for details. 

The projector can be chosen at will. If it is taken to be aj^^^a"'* it selects the subset of states (present in the Unruh) 
which contains that asymptotic quantum without specifying when it is detected. One can also consider the other 
limit in which one only specifies the moment of detection uq- In that case, Eq. (|25p is identical to Eq. (1^^ . up to 
an overal constant factor, as can be seen by taking /(r) = 6{t - uq) in Eqs.(55, 56) of jl^. From this we learn the 
physical meaning of Eq. ([24]): it gives the mean value of Trr when the in state is Unruh vacuum, and when a particle 
is detected on 1+ at u = uq. 

When (r, v) is also onX^ and coordinated by u, Eq. ([M)) is given by Eq. ([2]) (times 4-^), as can be seen using Eq. (|21|) . 
and Ku = —Ih^—kUk)- In this we recover the thermal correlations j3|] of the asymptotic radiation. When (r, w) is 
taken on the opposite null infinity , or sufficiently far away from the horizon so that w is constant, frr behaves as 
Eq. when using the mirror coordinate ku = IhIkUk)- Moreover, when (r, u) is near the horizon, K{r — rh) <C 1, 
and parameterized hy Uk, Trr behaves as in Eq. ([T]), as can be seen using Eq. (IT5|) and dllK/dr ^ — 2 for r = r/j. In 
this we recover that when probed near the horizon, for kx <^ 1, Unruh vacuum behaves like Minkowski vacuum. 

In addition to these three asymptotic behaviors, the non-trivial information contained in Trr of Eq. (j24p is the 
smooth interpolation from one to the other, which is represented in Fig. 1. From this we clearly see the gradual 
emergence from v = —8 of the energy flux associated with the partner on the other side of the horizon. What is non- 
trivial is the following. In the past of that time, Trr is essentially constant along the outgoing null lines Uk = const. 
and, as could have been expected, behaves exactly as Tiju would in Minkowski. On the contrary, near the horizon 
and in the future, Trr behaves very differently since the lines Trr = const, cross the horizon. This peculiar behavior 
could not have be found had we studied the 2pt function of Eq. ([7]) because the latter obeys Eq. ((Ij and depends only 
on U, even in a black hole geometry. This establishes that the use of Trr of Eq. (IT2]) with r affine is truly necessary. It 
should be also noticed that the above behavior of Trr cannot be found in Minkowski using affine coordinates either. In 
fact, it is characteristic of pair creation processes, as can be seen by comparing Fig.l to Fig. 1.1 of jll which describes 
pair creation in an electric field. As noticed in |10l | , th ese properties of Trr provide a clear answer to the long standing 
question: where is a Hawking quantum "born" ? |25| . 



2. Equal time correlations 



We consider Eq. ([18]) at equal EF time v = vq. Since Uk = e 



'Uk, the correlator 




(26) 
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Figure 2: Equal time correlations. On the left, we represent Eq. (|26|l. and on the right, Eq. (|29|l . both in the metric of Fig. 1. 
The horizon is at x,xo = 0. On the left, the signal diverges for x ^ xo whereas the pattern along x + xq = represents 
Eq. (|28|l . It is translation invariant in a; — a;o once w has reached a constant. On the right, the subtracted Sk is everywhere 
finite and regular. It is dominated by the correlations across the horizon. The subdominant patterns centered along x — xo = 
are due to the fact that the un-subtracted correlator in Eq. p9p decreases faster than the subtraction. On the diagonal, 
Sk{x,x) = {Trrix))^'^, the renormalized flux of Eq. ([T6)) . 



is given by Eq. (1^ with Uk replaced by Uk- Using Uk of Eq. (1^ . C/f (r, tq) diverges as ^ (r — tq) when r — tq, 
as expected since the Unruh vacuum is a Hadamard state. Moreover, when r — > oo, since Uk ^ — e^^"'', one has 

^^('^''^°)= 1. 2 ■ ^- (27) 

lOTT^ smh {K,{r — rQ)) 

The thermal noise of Eq. ([2]) associated with the asymptotic radiation is properly encoded in Ck since at fixed v, 
du = —2dr. Similarly, when x and a:o are on opposite sides of the horizon, since Uk ^ — e"^**!^! for x — ^ —oo, one 
asymptotically gets 

CK{r,ro)^ ^ -, (28) 

lOTT COSn [KX + KXo)) 

thereby recovering Eq. ([3]), and making contact with p^.[l9l|. see Fig. 2 on the left. To further investigate the gradual 
change from near horizon configurations to long distance correlations, we consider the subtracted 2pt function'* 



SK{r,ro) = drdr, 



{cf>{r)4,{ro)}K~^Hr-ro) 
in 



(29) 



The substraction term is the same as in Eq. ((T7)) . but considered here for all values of r — rg. From the the right 
plot in Fig. 2, we see that the residual signal is free of UV divergences and describes the emergence of the pairs on 
distances K(r — rh) = Ka; ~ 1, i.e. characterized by the geometry. 

One first notices that the correlations across the horizon in Sk are negative (as in inflation [1^). In fact the 
correlator {dr(j) dra(i>) k is negative "everywhere" , as can be seen from Eq. ((T9)) and Eq. (|22)l . We have added quotation 
marks because this correlator is a distribution. In fact, the coincidence point limit is ruled by the ze, see Eq. ([1]). It 
specifies that (the real part of) {dr(t> drg(j)) k diverges positively for r — > tq same v, and it ensures that the integral 
dU{du(j)d{jg(j)) vanishes in the Minkowski vacuum. (This is reminiscent of the behavior of {Tuu)fr., the mean 
flux evaluated in the Fulling- Rindler vacuum [1^.) We also notice that for close points, Sk is positive because 
the subtraction is larger than the "bare" term. Moreover when evaluated at the same point, Sk{x,x) is equal to 
{Trr{x)Y^^ of Eq. For a; ^ oo one verifies that Sk{x, x) = (—1/3) x Sk{x, —x). Finally we notice that Fig. 2. is 

symmetrical under x — > —x. This follows from the symmetry of 1 — in Eq. (|23p when k — R. When R/k ^ 1, this 
is lost, and 1 — resembles more to Schwarzschild. We hope to report on this case soon. 



* For simplicity we worked with the square root of Eq. I I26II rather than Ck itself. The reason is that the subtraction needed to obtain a 
finite expression for r — > ro is more complicated, as it requires three terms. The subtracted correlator of (TrrTrr) k possesses the same 
behavior as Sfc- 
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F. Time-dependent growth of flux and correlations 



So far we considered the stationary correlation patterns found in the Unruh vacuum. We now consider the early 
transient effects. As already said in Sec. III.B, they are not universal. However their late time behavior is when the 
state contains no high frequency excitations, i.e. when it is an Hadamard state. There is a simple and efficient way to 
characterize this behavior. It consists in assuming that the initial state specified at v = is the Minkowski vacuum. 
(This state can be found in a gravitational collapse when the infalling matter is a light-like thin shell [lOl-) In this 
state, V = fin, the 2pt function at is G""(r, ^o) = — + ie), as in Eq. ()17|) . In the future, one has 

G'"(r, v; ro,vo) = --^ ln(X'"(r, v) ~ X'"(ro, vo) + te), (30) 
47r 

where X™{r, v) gives the value of x = r — hit by the outgoing null ray issued from r, v when it crosses Vin- In this 
state, the mean flux of Eq. ([T6| becomes 



and the correlator of Eq. (|T8|) is 

(T„(.,r)T„(.o,ro))i„=(^^) {^^ (^X^.^,^,) - x\vo,ro) + ^e)^) ' ^^'^ 

In both expressions one has simply replaced Uk, the Kruskal coordinate encoding the Unruh vacuum, by X'" which 
encodes the time-dependent state which is vacuum at Vm- 

In the metrics of Eq. ([25)1 . these expressions flow towards Eq. (ITC]) and Eq. p^ . exponentially in v — Vin in the near 
horizon region where drW = k, and linearly when drW ~ 0. To show this we make use of hlK{x), solution of Eq. (|20|) . 
Calling X{Uk) the inverse function, and using Uk{x, v) = UK{x)e'^'", we obtain 

X"(r, v)=X \Uk{x) e-^^"-"-)] , (33) 

for all profiles w'^{x). To give a simple example, we use Eq. ()22p with R = k, and we get 

2«;X'"(r, v) = arcsinh [e.-"^'"-'''^'^^ sinh(2Ka;)) , 
dX" cosh(2Kx) 



dx [e2«(---.n) +sinh2(2K2;)]'/' 



(34) 

In Fig. 3, on the left, we represent Eq. ([32]) evaluated at equal time for different values of v — v-m- In these plots, 
one clearly sees the growth of the correlations across the horizon with a rate given by k [l^. One also observes a 
narrowing of the spread of the dominant correlations centered along x — xq. This is due to the progressive replacement 
of vacuum correlations ~ l/(x — xqY by the thermal ones '-^ / sinh^ k,{x — xq) of Eq. (P7|) . On the right plots, the 
subtracted function 5'in(r, tq), the equivalent of Eq. (p9)) evaluated in the in vacuum, displays both the growth of the 
long distance correlations, and the modifications of the local correlations. 

It is worth analyzing these time dependent effects through two other perspectives. In Fig. 4, on the left, we present 
the mean flux of Eq. pip in the x^v plane from the onset of the vacuum at f = fin ~ 0. {Trr)in contains transients 
which propagate along u = w — 2a: ~ 0. It then reaches a constant x-dependent proflle. For a; — > cxi one recovers the 
standard value K^/127r [9]. In the metric of Eq. (1^^ with k = k, {Trr)in crosses for sinh2a; = ^/2, i.e. x ~ ±.57, 
and on the horizon, it is negative and equal to — 2K^/127r. When including gravitational back-reaction effects, this 
term participates to the evaporation of the black hole, see [l^ for a numerical analysis in similar settings. On the 
right plot, we show fuu{u)\uo of Eq. (p4)) . the flux correlated to a late detection at uq on Xj^, that we evaluate on 
and parameterize with u, the mirror coordinate ku = IukUk- There is no correlations from to for uq < — 1 
in conformity to the fact that no Hawking radiation as yet reached the null infinites. Then for positive uq, the 
correlations settle to a stationary pattern centered around Uq + u = found in the Unruh vacuum. 

In conclusion it is interesting to observe that, even though the transients give rise to a higher value of the mean flux, 
they are not associated with stronger correlations across the horizon, and this because, unlike the steady Hawking 
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Figure 3: The growth of equal time correlations. We represent Eq. (|32p at equal time on the left, and Eq. (|29p on the right, 
after a lapse of time «(w — Wn) = 1 in the upper plots, and k,(v — win) = 4 in the lower ones. On the left, one observes the 
growth of the correlations across the horizon centered along x + xq =0. One also observes a narrowing of the correlations 
centered along x = xq. On the right, Si-n{r,ra) displays two distinct features. A strong signal associated with the building up 
of the correlations across the horizon, and subdominant patterns on both sides of the horizon due to the growth of thermal 
correlations of Eq. (|27|) . At late times, both patterns asymptote to the stationary ones of Fig. 2. 




1 2 3 4 -1 0123456 



Figure 4: Time dependence of the mean flux and of the asymptotic correlations. On the left, we represent Eq. (|31|l with Eq. (|34p 
from the onset of the vacuum at w = till nv = 8, and for k,x from to 4 (It is symmetric under x — > —x). The transients 
propagate on null lines u = v — 2x ^ Q. After they passed, {Trr)in is u-independent. On the right plot, we represent Tuuiu)\uo 
of Eq. (|24|) in the uo, m plane with u vertical, with uo and U defined on the future null infinity vq = v = oo. For uq < 0, before 
the transients, there is no correlations across the horizon. Instead for uo > 3, Tuu{u)\uq only depends of wo + it, and is given 
by Eq. (O. 
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radiation, the transients are not composed of entangled pairs of opposite frequency w. This can be checked by 
comparing Ty^u{u)\u to the asymptotic flux {Tuu)in evaluated for x,?; — >■ oo, with u = v — 2x fixed: 

(7Mti(u))i,j = ToZT-, T 7TT2' 
487r (1 4- g-2Kti uy 

T^n(S)U„ = f^j ((e2-'° + l/4)(e-2- + i/4)(arcsinh(e-™V2)+arcsinh(e+'=72))') . (35) 



IV. DISPERSIVE THEORIES AND ANALOGUE BLACK HOLES 



In a non-homogeneous medium, linear density fluctuations obey a relativistic d'Alembert equation in a curved 
metric when their wave lengths are larger than the inter-atomic distance [l3|- Instead, for shorter wave lengths, the 
propagation becomes dispersive Assuming the speed of sound is constant and set to 1, it can be described by 

n'^ = {uj-wpf =F'^{p^), (36) 

where 17 is the co-moving frequency measured w.r.t. to the atoms, w the velocity of the fluid, and p, uj the wave vector 
and the frequency measured in the lab. The dispersion can be incorporated in a modified field equation [ll| 

{dr + d^w) {dr + wd^) - ± = o, (37) 

which reduces to the relativistic one when sending the dispersive scale A — )■ oo. For simplicity, we chose quartic, super 
(-f-) or subluminal (-) dispersions: — ±p'*/A^. 

Eq. (j37p can then be used to study the impact of dispersion on Hawking radiation. In [ll| . the thermicity and 
the stationarity of the asymptotic radiation have been shown to be robust, i.e. hardly affected by dispersion when 
K <C A. This is sufficient for recovering Eq. (j27|) . By constructing wave packets of in modes, it was then shown fl^ 
that at large distance from the horizon one also recovers the correlations between Hawking quanta and their partners, 
and this, even though the early propagation was radically affected by dispersion. This second aspect is sufficient for 
obtaining Eq. ([28l) . Fig. 2, and the late time properties of Fig. 1. 

To clarify this, we shall compare the properties of wave packets of 0™, the in modes of Eq. ([57]). 

^^(t,x) = J du;e-'^^cb:;{x)U, (38) 

where fui selects the wave packet, with the correlation function 

G™(r, a;;ro,a;o) = ((/)(t, x) (/)(to, xq))™ 

d..e-^"^C(^)[e"^"^°CK)]*, (39) 



evaluated in the stationary in vacuum. When dispersion is weak, the similarity of expressions guarantees that 
similar patterns will be found. However, since dispersion grows as approaching the horizon, the way one probes 
the correlations, i.e. by extracting some limited range of lo through f^j in Eq. ([38|) . or not as in Eq. (f39|), can lead 
to different behaviors. We now review the relevant points for achieving this comparison following [la . |30|; other 
treatments are mentioned in 13111. 



A. Kinematics 

From a relativistic point of view, the presence of dispersion defines a preferred frame fl3| which allows to define 
new scalars. This is best seen by " covariantizing" Eq. ([36]), i.e. by introducing a unit time- like vector field u^, and 
viewing the field (j) as propagating on a manifold endowed with both the metric and u^. Then the energy in the 
preferred frame, and the spatial momentum perpendicular to it are respectively 

n^u'^p^, p^s^p^, (40) 

where —u^ = 1 = and s^u^^ = 0. The PG coordinates used in Eq. (1571) can then be invariantly defined by dx = s^d^ 
and c?T = Uf^dx". The field u'^ also defines the scalar density p = u^u'^T^j^n which corresponds to the proper energy 
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that observers following dx^ /dr = would measure. It is interesting to note that in the hydrodynamical limit one 
gets 

P = T,x^ {d^4>?- (41) 

In other words, p coincides with Trr of Eq. ([T2l) with r defined by Eq. (fTTj) . This non-trivial correspondence follows 
from the affinity of r at fixed v and r, see the remarks after Eq. ([SJ. 

In the same spirit, we notice that the surface gravity measured with respect to the preferred frame is also scalar. 
It is given by the expansion (32j = u^^ — dxW evaluated on the horizon. (The second expression is valid in PG 
coordinates). The ambiguity of the scale of the surface gravity in covariant theories is thus removed when using . 



B. The modified modes 



In linear field theories, the modifications of the stress energy correlations due to dispersion will stem from the 
modifications of the modes, solutions of Eq. (1571) . When A ^ k, these are modifications localized near the horizon, 
for K(r — r/i) = <C 1. As noticed in [l^ it is appropriate to work in the p-representation (with x = idp) with w{x) 
linearized: w = — \ + kx = —1 + indp. In this representation, Eq. (|37|) becomes 

{uj - pw) {uj - wp) (t)uj = (f)^. (42) 
Using w — —1 + indp, the modified modes have the form, for details see the Appendix of [30| . 

^^^4>l^ e-'^/^x{p), (43) 
where — is the standard dispersion less mode, and where x obeys 

-n^d^x^^X^H^X- (44) 

The in state which generalizes the notion of the Unruh vacuum, the " free falling vacuum" , is characterized by the 
positive norm modes which contain only positive 51 of Eq. PO]) . These in modes are related to the Unruh modes, see 
Eq. (El, 

by M 

C = X e-^P^^x{p), (46) 

where x{p) is the solution of Eq. (pi)) with a Wronskian equal to x*dpX ~ X^pX* = 'ii/n. The corresponding WKB 
solution is 

xiP) = JJ^ cxp (^z 1^ H{p')dp'/n^ . (47) 

It provides a reliable approximation when A ^ k. In the limit A — > oo fixed p, — > 1 and e^'^/'^x — > 1, thereby 
implying that 0™ smoothly gives back (p^.^ 

To understand the impact of dispersion, we now study the characteristics of Eq. (1371) since the maximum of 
correlations will be localized along them. Having already the modes in p-space, the simplest way to get them is to 
consider Eq. (|38p in p-space, and look for the stationary phase condition in w. Using Eq. (|46|) . since x is independent 
of 07, one gets 

p(T)=poe-'^^ (48) 



^ At fixed X instead, tlie limit A — > oo can be singular as some roots p{x,u}) of Eq. I|36p are sent to infinity. The p-WKB approximation 
sfiould not be confused with the usual one defined in x-space. For the Airy function, the modes in p are exactly given by their p-WKB 
approximation. Similarly here, the corrections to Eq. I|47| l are negligible when A/k 2> 1. This has been confirmed by numerical analysis, 
see [35ll for a detailed study. From now on we neglect them and work in the adiabatic approximation with Eq. I I47I I. 
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irrespectively of the dispersion relation F, and thus as in relativistic theories. 

To get the modified characteristics in x, we use uj — wp — F, the root of Eq. p6p describing the right moving modes. 
Using w ^ —1 + KX, one gets (s^l 

2 p 

'^[xUp)+x-Up)] = l-H{p). (49) 

The first equation is again independent of F and coincides what is found in relativistic theories, namely, when 
propagated backwards in time, pairs of characteristics pill up exponentially in PG time. The second equation tells us 
that the "center of mass" of a pair which is centered on the relativistic horizon for 1 — ^ 1, i.e. p <^ A, gradually 
moves away as p increases. For sub (super) luminal dispersion, H < 1 {H > 1), the pair is sent outwards (inwards). 
For quartic dispersion the momentum at the turning point is given by ^ = 2A^|a;| [34i] . Because of this movement 
away from the horizon, the increase of p and the focusing of x^, — x-^ will stop when the pair reaches \kx\ ~ D, where 
the gradient dxW drops down. A straightforward calculation gives that the focusing stops for p ~ KD^/^. 

The end of the focusing and the movement away from the horizon are the principal consequences of dispersion. 
They imply that the early properties of Fig. 1 are inevitably modified, as we shall see below. 



C. Correlations from wave packets 
To show how dispersion affects the correlation pattern encoded in Eq. ([38]) . we need Eq. ((46|) in x-space 

Far away from the turning point, one can evaluate this integral at the saddle point approximation (since it is reliable [igI . 
[s^p. and decompose Eq. (|50p in terms of outgoing modes defined for low momenta. Doing so one finds 

C {x) = [e{x) {x) + e{~x) [ip.^ (x))*] , (51) 

where ipuj{x) is the x-WKB mode of Eq. (|37p with unit norm. Explicitely it is given by 



t^-. _ , ^pUx) exp (» /;„ dx'pUx')) 

where Puj{x) is the corresponding low momentum root of Eq. p6l) . One easily verifies that these out modes are 
identical to the relativistic ones for p^^ ^ A. Moreover, for k <C A, up to a phase, one finds z^^ = e"'^"/". Hence 
Eq. ([51]) gives the equivalent of Eq. (K2\i . ^ 

In addition to the above two low momentum modes, there is a third saddle -on the left (right) of the horizon for 
super (sub) luminal dispersion, in conformity with Eq. (j49p - which gives a high momentum mode. Its WKB wave 
is also given by Eq. ([52]) with p^j being the unique large positive real root of Eq. ([36)) . One verifies that its overall 
coefficient is unity in conformity with the fact that it describes the incoming mode that shall be scattered. 

Therefore, considering Eq. psp with fi^ centered around < (D ^ Wmax, we get two results. First, at late times, 
using Eq. ((ST|) and Eq. (|52l) . one finds two low momentum packets following Eq. with uj — ±a), where the negative 
frequency packet has its amplitude reduced by Za, as in Eq. (IA3|) . Since z^ = e~^'^/", and since both ipi^ of Eq. ((5T|) 
behave as relativistic modes once p is small enough (F — p <ti p), at large distances, the pattern is indistinguishable 
from the relativistic one obtained by replacing Eq. (|46l) by Eq. psj) . 

Second, at early times, only the incoming high momentum mode constructively interferes. It has a mean positive 
frequency w, follows the second line of Eq. (1491) with p ^ Q, and leaves the near horizon region with p A. This is 
completely different from what is obtained in the relativistic case. Indeed, using Eq. (1451) , p would keep increasing for 
ever following Eq. (j48|) . and the spread in x correspondingly decrease as 1/p. 



This approximation is valid provided oj is sufficiently small. For quartic dispersion there is a critical frequency ojniax, related to both 
A and the asymptotic velocities w{±od), above which identically vanishes 35]. From the numerical results of that ref., a good fit is 
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From the analysis of wave packets, we have thus reached two important resuhs. On one hand, the low momentum 
(late time) properties of the relativistic pattern of Fig. 1 is unaffected by dispersion. On the other hand, the early 
properties of this pattern will be radically affected by dispersion since the peak of correlations will follow Eq. (|49|) . 
as represented in Fig. 4 of Tg']. To further investigate how dispersion affects the properties of Figs. 1 and 2, we now 
consider the pattern encoded in Eq. (|39)) rather than in Eq. ([38)) . 



D. Correlations in energy density 

We start with the correlations of a relativistic field expressed in the present language. Since fl = p, the correlation 
function of p of Eq. (|4T|) is 



{pix,T)pixo,To)),r, = [dA,p dc^e-^"(^-^")G™(x,^o)) , 



(53) 



where G]^ = </)^^(x)(0^"'(xo))* is the uj component of G^^ of Eq. ffSO]) . In the ^-representation, using Eq. f|45|l . one gets 



oo 

iujT r^in ( 



J — OO 

- e{p)e{pa)-^5{KT + \n(p/po)). (54) 

On one hand, we recover the classical evolution law of Eq. P5|) . On the other we learn that in the Unruh- vacuum, 
at r = To, only configurations with equal values of p contribute. There is no spread in p in this state. When inverse 
Fourier transform, one gets 

Gk{x, t; xo,Q) ~ — ln(a; — xq e""^ + ie). (55) 
An 



We recover the standard result, Eq. (IB5|) . expressed in PG coordinates. At equal PG time, we notice also that the 
argument of the log is x — xq. We shall return to this point below. When computing dxdxgGK one obtains 

1 e'^^ 

dxdxoGKix,T;xo,0) = - — - , (56) 

in [x — Xq e"^ + ze)^ 



which is the square root of Eq. (|T9)) in the near horizon region where Uk ^ ~x. 

When introducing dispersion, p receives corrections with respect to {dxipY due to the non-linearities of . This is 
hardly relevant for us, because in the near horizon region, the momenta p are much smaller than the UV scale A. They 
are of course modifications in the UV sector of the theory, but these ultra local effects are the same as in Minkowski 
space. Therefore the main modifications will come from the replacement of the Unruh modes by the modified ones. 
Using Eq. (gH]), Eq. dM]) is replaced by 

expi P [i?(y) - lldp'/K 
G»(,,r;po,0) = G.(,,r;,o,0) x - (57) 

It should be noticed that the frequency Wmax mentioned in footnote |6l will cut out the integral in Eq. ([53]) , below 
^i^max for superluminal dispersion, and above w^ax for sub luminal. We ignored for this UV cutoff in computing 
Eq. (|57p because Eq. (|46p is no longer trustworthy anyway when uj — > Wmax- In fact, in the adiabatic approximation 
of Eq. (|47p. there are cancelling errors, in that the next equation can be shown to be exact. 



1. Equal time correlations 

At equal times, Eq. ([57| gives 



G™(p,po,^r = 0) = (58) 



We proceed as in [Tsll . Nevertheless the forthcoming equations differ in several respects. 



16 



and in the a;-representation, one has 



G»(.,.o,^r.O) ^ / (59) 



The only effect of dispersion is to replace in the denominator the relativistic law fl — p hy ^ = F(p). Therefore 
Eq. (|58p is exactly what one obtains in Minkowski vacuum in the preferred frame. The reason is again that x is afHne 
at fixed t: ds^ — dx^ . In fact, as in the relativistic case, see Sec. III.F, a non-stationary vacuum can be defined at a 
given time but for all values of x by plane waves e'^^ with p > 0. Then, the negligible character of the non-adiabatic 
corrections to Eq. (j47p in the near horizon region gives Eq. (j57p which implies that in that region but at all times the 
vacuum stays characterized by p > 0.® This also implies that the dispersive version of the time-dependent Eq. (|30|) 
will evolve towards the stationary Eq. ([5^ . as Eq. ([50]) evolved into the stationary function in the Unruh vacuum. 

From Eq. ([5^ several consequences can be drawn. If one probes the in state for k\x\ <^ 1, the deviations w.r.t. to 
the relativistic case ioi x — xq < 1 /A {p > A) are the same as in Minkowski, and are therefore insensitive to presence of 
the black hole. If one probes the in state further away from the horizon and for momenta p < A, since ^ p, Eq. (j39p 
will behave as the relativistic function, as it obeys the same equation, and possesses the same initial conditions. Hence 
the whole analysis of Sec. III.E.2 applies. In particular, as soon as w is constant, dxW ^ k, Eq. ([55)) will obey Eq. (P7)) 
when both points are on the same side of the horizon, and Eq. (j28p when one is on either side. ® Thus the properties 
of Fig. 2 (left) are not affected by dispersion when k ^ Wmax- Those of Fig. 2 (right) are not either when subtracting 
the dispersive expression that replaces the log in Eq. (P^. because Sk varies on scales 1/k ^ 1/A. 

The insensitivity of Fig. 2 against introducing dispersion is quite surprising since, as discussed before, we expect that 
the properties of Fig. 1 be affected by the drift of Eq. p9)) which occurs for rather low momenta ^ A^/^k^/^ ^ A. 
The reason of the disappearance of the drift (at equal PG time) is the following. In Eq. (|55)) . because we are summing 
over w, we erase the coherence in x-space that exists in each uj sector, thereby recovering the translation invariance of 
the in state, as in Eq. (jB5l) . In other words, it is only when isolating some uj content out of all vacuum configurations 
that the early pattern characteristic of wave packets emerges. This deserves further comments. 

Given Eq. (|59p , what can be said about the entanglement entropy ? The regular behavior of the dispersive in modes 
and the entanglement in Fock space between states of opposite ui, see Eq. (jB2l) . were exploited in [30] to argue that 
the entanglement entropy of a black hole is finite (in 1+1 dimensions). However, using Eq. (|59p one would conclude 
that upon tracing over inside configurations a; < 0, one would obtain the same (diverging) result as in Minkowski [37| . 
This conflictual result indicates that there is probably no unique notion of the entanglement entropy. Therefore to 
get a well defined result, it is needed to specify what one exactly means by "tracing over the inside configurations". 

We saw that the 2pt correlation function (at equal time) does not display the characteristic pattern of wave packets 
with a given frequency content. This is quite general. It was discussed in J6| when studying the correlations amongst 
particles emitted by accelerator mirrors, and in a inflationary context in [26 

2. Correlations at different times 
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When T ^ 0, in the relativistic case one gets Eq. ([55)1 . Instead, Eq. ([57| gives 

where S — x — xqc'^'^ . The non-trivial modifications of the correlations due to dispersion are best seen by evaluating 
this integral at the saddle point approximation. The value of the saddle answers the classical question: given that 
one starts at a;o; '''o = smd ends at a;, r, what is the momentum at that time ? It is given by 

- [e«^(if(p,e^^) - 1) - (HM - 1)] . (61) 



Eq. I I58I I is also obtained in the adiabatic approximation, in cosmological backgrounds when the preferred frame is aUgned along the 
cosmic frame. Moreover, this correspondence becomes exact (beyond the adiabatic approximation) when considering de Sitter space when 
the Hubble parameter H = k since the linearized expression w = —1 + kx describes this space in PG coordinates when — oo < x < oo. 
In this we recover what has been found in Bose Einstein condensates when looking at the density-density correlation function 
In that case, in the hydrodynamical limit, the atom density fluctuation is given by dx(p, and the correlation corresponds to Eq. I I56II . 
These remarks raise the question of the choice of (the set of) observables used to probe a quantum state. To give a concrete ex.: in 
inflationary cosmology, it is generally assumed that the large amplification experienced by primordial fluctuations erase all quantum 
properties and would give a state indistinguishable from a stochastic ensemble of classical fluctuations. In [38l | it was shown that 
irrespectively of the amplification there exist observables exhibiting violations of Bell inequalities (for linearized modes). 
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To understand the implications of this expression, we consider three regimes. First, if K(5 <C 1 and kt ^ 1, the two 
points are ahxiost connected by a null ray and the red-shifting effect is moderate. In this case J3,/A <^ 1 and one 
recovers the relativistic behavior of Eq. (f55|) . 

Second, we study the non-trivial correlation far away from the light cone (but still in the near horizon region). For 
definiteness we restrict attention to quartic laws = ±p^/A^. We expand Eq. ((6T|) to first order in 1/A^, and 
using H — 1 ^ ±p^/2A^, we get 

^<5 = ±1^(6^-^-1) (62) 

When a; = Xq > 0, irrespectively of the sign of r, there is no (real) saddle for the -I- sign, i.e., superluminal dispersion, 
in agreement with Eq. (|49p which says that both partners are dragged inside the black hole horizon. Instead, for 
subluminal dispersion, since they are both dragged outside, there must be a non trivial solution. To confirm this, we 
take KT such that e'*'^ ^ e^'^'^ . In this regime, Eq. (|62p reduces to koiq = ple^'^'^ /2h? . From this we can deduce a;, 
the mean value of the frequency corresponding to the trajectory that goes from xq back to it in a lapse equal to r. It 
is approximatively given by 

L^^^y/2K{KXQfl^ e"^^ . (63) 

This result can also be derived using Eq. (H^ (and applied to superluminal dispersion for x < 0). Thus, when studying 
G{x,T;X(j,Q) at sufficiently large kt, unlike what we found in Eq. (j59p . the correlations are now in agreement with 
the locus of constructive interferences of wave packets because only a limited range of frequencies centered about 
significantly contributes. This confirms that near horizon behavior of Eq. (|53p will completely differ from that of Fig. 
1, and will be similar to those of Fig. 1. of [30|. What remains to be clarified concerns the profile of Eq. ([55)1 at early 
times. Namely, at fixed a;o,To, what is the trajectory of the maximum of Eq. (j53p . and what is its spread in a; as a 
function of xq, r — tq, and A ? We conjecture that both of these quantities are ruled by Wmax of footnote [51 

What can be studied |3l is the " off-shell" limit of very large blue-shift e'*'^ 3> 1 encoded in a backward propagation 
at fixed X and fixed 5 = x — xqc'^'^ . This limit displays how dispersion tames the " trans-Planckian" behavior found 
for the relativistic field. In that case, Eq. ([55)) gives 

1 e'^'^ 

dxd,j:^,GKix,0]Xo,-T) ^ - — (64) 

An [0 + ley 

how ever large is kt > 0, in agreement with Fig. 1. In the dispersive case, using Eq. ([5U1) and H{p) — 1 ^^.^ ±p^/2A^, 
one has 

When the blue shift is moderate, i.e. e'^'^ < k6 (A/k)^/"^, Eq. ([55]) behaves as Eq. ([M)) plus corrections in e^'^'^ /K^kS^ <^ 
1 that can be computed perturbatively, as can be seen by changing variable p ^ q = pS. Instead, when the blue shift 
is large: e'^^ > kS (A/k)^/'^, the integral becomes independent of S as is seen by using k — pe 

''^/(A2k)1/3. Explicitly 

one finds 

9,9,„G^"(a;,0;xo,-r) ^ (A^Kf/^e-- x C±{k/A,t), (66) 

where C±(k/A, r) are slowly varying functions which stay bounded for t oo. This exponentially decreasing result 
can be seen as the contribution on the horizon of the tail of the configurations with high p which follow the second 
equation in Eq. 

This smoothing out of the relativistic behavior is very reminiscent to what was found in [21] when studying the 
backwards evolution of G™{x,0;xo,—t) of a relativistic field propagating in a stochastically fiuctuating black hole 
metric. In addition, for nearby points, G™ in a stochastic geometry also behaved as Eq. ([5^ . as can be seen in Eq. 
(4.7). Based on this similarity it was argued j23 [ that when taking into account the gravitational radiative corrections, 
the dressed Green functions should effectively behave near a black hole horizon as in Eq. ([66)) . thereby reinforcing the 
idea that the unbounded growth of Eq. ([Ml) cannot " accommodate gravitational non-linearities" . 

V. CONCLUSIONS 

We showed that the monotonic energy correlations found in the vacuum Eq. ([T]) gives rise to a maximum of 
correlation across a Rindler horizon when re-expressed in terms of coordinates associated with accelerated systems. 
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see Eq. ([3]). This maximum is not a mere coordinate artefact as it affects the combined state of co-accelerating 
systems. 

In Sec. Ill, we transposed this analysis to stationary black hole geometries, and recalled that the regularity of the 
state across the horizon and the inertial character of asymptotic observers are essential to provide a physical meaning 
to the thermal correlations of Eq. ([2]). When considering black hole geometries which contain asymptotic regions on 
both sides of the horizons, the correlations of Eq. ([3]) are found at large distances when using inertial coordinates. We 
then make use of the affine parameter of Eq. (jlip to obtain an invariant description of the energy correlations in 
the entire space-time. We compared the correlation pattern associated with a late detection. Fig. 1, to that obtained 
at equal EF time. Fig. 2. In both cases the gradual emergence of a maximum of correlations across the horizon is 
clearly visible. By considering the subtracted correlations of Eq. (1^^ . we saw that the remaining signal is dominated 
by the long distance correlations across the horizon, and also contains a sub-dominant local contribution associated 
with Eq. This analysis was generalized in III.F. by including the transients effects which precede the stationary 
patterns found in the Unruh vacuum. 

In Sec. IV. by studying both wave packets and correlation functions, we studied how these patterns are modified 
by dispersion. Far away from the horizon, the pattern is robust, i.e. hardly affected by dispersion. Close to the 
horizon we saw that dispersive effects show up differently depending on how one probes the state. When probed at 
equal PG time, the correlation function is translation invariant, and as in Minkowski, see Eq. (|59|) . Instead wave 
packets of in modes centered around a given frequency uj display a characteristic pattern which follows the modified 
characteristics of Eq. (1491) . When the momentum has sufficiently increased (in a backward in time propagation) the 
wave packets are dragged away from the horizon, and, the blue shift effect saturates. This behavior is recovered from 
the correlation function when considered at different times and appears through a non-trivial saddle point in Eq. (j60l) . 
When considering the correlation function for two points separated by a very large PG time, the drag w.r.t. the 
relativistic horizon results in an exponentially suppressed amplitude in the place of the exponentially growing result 
found in relativistic theories, compare Eq. (|64p with Eq. (j65|) . These properties are reminiscent to what was found 
when considering field propagation in a fluctuating black hole metric, and could possibly be found when taking into 
account gravitational interactions at the quantum level. 
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Appendix A: Unruh modes 

Firstly, they are solutions of d'Alembert equation dudv4>ui = 0, and thus only depend on either U or V . Secondly, 
they have a fixed boost frequency uj, i.e. they are eigenmodes of 

idu(l>uj = -iaUdu4>uj = ^'t'u- (Al) 

Thirdly they are only composed of the positive norm modes: (j)^ — e^*^'^/(47r51)^/^ with > 0. Explicity they are 
given by 

6l(-C/)(-at/)^"/" + X e{U){aUy'^'''] , (A2) 



(47rw)i/2 

where the normalization obeys |a(^P = (1 — e-STrt^/a-j-i^ ^^^^ where z^j = e~'^"/°. This factor arises from the ie 
prescription which specifies that the analytic continuation from the R to L quadrant must be done in lower half 
complex U plane. As in Eq. (jT]) this prescription comes from the fact that only positive frequency fl = idu modes 
contribute. 

Moreover, they are globally defined, —oo < U < oo, and form a complete and orthonormal basis of positive norm 
modes when — oo < u < oo (with respect to the standard Klein-Gordon product). Hence the Minkowski vacuum can 
be alternatively defined as the state annihilated by the destruction operators Quj associated with these modes. 

Thus, when a quantum system is (linearly) coupled to (j) which is initially in the vacuum, the transition amplitudes 
will contain some (linear) combination of the (j)^. When the system is not accelerated (e.g. inertial), the decomposition 
[2|) presents no interest since the system will cross U = 0. On the contrary, when it is uniformly accelerated in, say. 
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the R quadrant, Eq. (jA2l) guarantees that every transition occurring in the Minkowski vacuum defines a partner wave 
in L, see App. C. for more details. 

This R — L partnership can be studied in simpler terms and without referring to accelerated systems by constructing 
wave packets of Unruh modes 

POO 

0= / dLjf^cl>^^^'' + ^'^. (A3) 
Jo 

Eq. ()A2|) thus implies that to every packet (j)^ localized in R will correspond 0^, its partner wave in L. More can 
be said: since z^i in the r.h.s. is real for all lu, when (p^ constructively interferes around some Ur < 0, Eq. (|A2p 
guarantees that (j)^ will do so near —Un. This explains why the maximum in Eq. arises for opposite values of U. 
In addition, from the fact that high uj are exponentially suppressed by = ^-'^^1°- ^ the maximum in Eq. (|3]) cannot 
diverge as it does in Eq. ([2|). Instead it must scale as a"*. In brief, the mathematical properties of Eq. ([3]) are deeply 
rooted to those of the modes i/);^. 

There exists an efficient way to encode the properties of the Unruh modes which turns out to be very useful when 
analyzing Hawking radiation in the presence of dispersion. It consists in computing the Fourier transform at fixed t: 
(f'uip) = J dz e~^P^(l)i^/{2ny/^ . Taking into account the ie in Eq. (jA2[) . up to an irrelevant phase, one finds 

The restriction to positive p follows from the fact that property only positive fl contribute to (j)^], and from the 
dispersion relation ft — p which describes right moving modes. We also note that when considering Eq. (jA3[) in the 
p-representation, the two wave packets found in the ^-representation on either side the horizon are now described by 
a single packet in p-space. This is characteristic of pair production phenomena, see e.g. Sec. 1.2-1.3 in j5|]. We shall 
return to this in Sec. IV. 



Appendix B: Fulling-Rindler states and partner-ship in Fock space 



If a mode analysis is sufficient to understand the behavior of Eq. Q in space-time, to have a deeper quantum 
mechanical understanding of Eq. ([3]) we analyze the bi-partite structure in the Fock space when using states with a 
fixed frequency ui. To this end, we should discuss yet another property of Eq. (|A2[) . It concerns the fact that the 
mode on the left {U > 0) has a negative norm (for oj > 0) thereby implying that the norm of the right component 
is correspondingly larger. This invites to consider the inequivalent quantization of (f) based on the "Fulling-Rindler" 
(FR) modes. These are normalized eigenmodes of frequency oj, Eq. (jAip . localyzed either on the right, or the left, of 
[7 = 0. Hence we re- write Eq. (|A2p as 



+ (3^ {(j}^y , forw>0. 



<l>u = a^u + P-u i^'lS ' for < 0, (Bl) 

where cj)^ — e^^^^ j (A.txlS)^^'^ {(j)^ — e^^'^^ / {AttloY/"^) vanishes on the left (right) of the horizon, and where jS^ — z^^auj- 
One easily verifies that — /^^ = 1, which implies = (e^'^"/" — 1)^^. For each w > 0, Eq. (|B1[) defines a (two-mode) 
Bogoliubov transformation relating (0^^, (t>-ui) to {(t>^, i/if^). This implies that the vacuum can be written as a product 
over w > of two-mode squeezed states 

|0) = H^ (— exp (z„ a^a^A \0)r |0)l, (B2) 



where the i?- vacuum \Q)r is the state annihilated by the afj, the destruction operators associated with the cf)^, and 
similarly for the L sector. Since the squeezing operator is quadratic and diagonal in a;, for free fields, all expectation 
values are expressible in terms of the following two VEV 

{{a^Ja^) = {(atVat) = W^f, (B3) 
(a^ai^) = l3^aZ^ z^\a^\\ (B4) 

It is now instructive to see how these two VEV enter in Eq. ^ and Eq. To this end, we consider the (U 
contribution of the) 2 point function of (f> 

{<P{U)(j)(Uo)) = I dujcj,^{U){MUo)y =-^HU ~Uo~ie). (B5) 
./-no 47r 
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Using Eq. (IB1[) two different expressions are obtained depending if both points are on one side, or on either side, of 
U = 0. Explicitely, when both U are negative and written as —all = e""" one has 

(0(C/) 0([/o)> = ^ (|a^|2e-''^("-"°) + |/3„|2e+^"("-"°)) . (B6) 

Instead when one point, say [/, is positive and written as all — e^^, one has 

(0(C/)0(;7o)> = ^2Re(|aj2^:e— (B7) 

One sees that Eq. (jB6p and Eq. ^ arise from diagonal terms, hence weighted by \[iuj\^ of Eq. (IB3p . whereas Eq. (jB7p 
and Eq. ([3]) arise from interfering terms weighted by a*/3i^ of Eq. (IB4[) which encodes the entanglement, in Fock space, 
amongst the R and L sectors. 

It is an interesting exercice to verify that when using the actual functions for a^j and /3w, Eq. (|B6P and Eq. (jB7P 
both give back, as they must, the log AU of Eq. (jBSp . Therefore, they are only complicated re-expressions of Eq. (jB5l) . 
Similarly, Eq. (jB2p is only a mathematical re-expression of the Minkowski vacuum. However, it prepares the analysis 
of the physical processes related to the Unruh effect, to the quantum fluxes emitted by non-uniform mirrors [3, 6|], 
and to black hole physics. In these three cases, there is an external agent -respectively an accelerated system, a 
non- uniform mirror, a non-trivial metric- which acts on the field and "transforms" the FR states into asymptotic 
states. This particularly neat in the case of the non-uniform mirror discussed in [3], see Eqs. (3.27-3.28), see also 
Section 2.5 in Q. 



Appendix C: The conditional value associated with a detection in 7? 

We recall how a detection of a FR quantum in R defines first, a partner state in L, and second, a projector which 
allows to define the conditional value of an operator associated with this detection. We describe the detected quantum 
in R by 

/•oo 

\^r) = / dwfUa!^)^\0)R. (CI) 



The EPR partner state is defined by reducing the bi-partite state. In the present case, the latter is the Minkowski 
vacuum expressed as Eq. (|B2I) . The partner state is thus 

POO 

- (*7?|0) ^ / dujz: fZ [a^mL. (C2) 



One notices that its Fourier components are fixed by 2;* and /*, i.e. by both the state and the complex conjugated 
of the component of selected wave packet. 

It is instructive to relate these two states to the (classical) wave packets of Eq. (jA3l) . To this end we introduce the 
projector IIr = |^_R)(5'ii:| and consider the value of (fP' conditional to the fact that the detection took place (for more 
details see Q) 

^2 = (O|02n^|O) ^(^2)fl^(^2)L_ (C3) 

A direct calculation gives 

/ POO \ / nOO 

= I / duJZzlc^^S- (C4) 

When = a^Juj where J^j given in Eq. (jA3p . (0^)^ exactly gives of that equation. Similarly the first factor 

in the first line is ^/j. The second factor is not its complex conjugated due to the presence of l^tjp in the integrand. 
However for wave packets with a small spread uj wrt a, this quantum mechanical feature (whose consequences are 
discussed in [To|) does not significantly affect the spatial properties of the R wave packet. Thus we basically recover 
the modulus square of the two packets of Eq. (|A3p . 



21 



The lesson of this exercise is that the pattern obtained by constructing wave packets of Unruh modes as in Eq. (jA3[) 
offers a rehable description of the quantum correlations across a Rindler horizon (when the spread in u is small enough) . 
This transposes in black hole metrics (without and with dispersion) and implies that the correlation patterns of in 
modes also offer a reliable description of the quantum correlations across the horizon. 

Finally, we mention that the study of highly excited coherent states, see App. C of [20], offers another way to 
relate the packets of Eq. (|A3p to quantum states. Using these coherent states, one can verify the agreement of both 
descriptions in describing the R — L correlations. 
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